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The magnetars are believed to be highly magnetized neutron stars having surface mag- 
netic field 10 14 — 10 15 G. It is believed that at the center, the magnetic field may be 
higher than that at the surface. We study the effect of the magnetic field on the neu- 
tron star matter. We model the nuclear matter with the relativistic mean field approach 
considering the possibility of appearance of hyperons at higher density. We find that the 
effect of magnetic field on the matter of neutron stars and hence on the mass-radius 
relation is important, when the central magnetic field is atleast of the order of 10 17 G. 
Very importantly, the effect of strong magnetic field reveals anisotropy to the system. 
Moreover, if the central field approaches 10 19 G, then the matter becomes unstable which 
limits the maximum magnetic field at the center of magnetars. 
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1. Introduction 

Anomalous X-ray pulsars and soft 7-ray repeaters are observationally identified with 
highly magnetized neutron stars, known as magnetars, with surface magnetic field 
- f0 14 - 10 15 G. m The processes of supernova collapse will leave behind a strongly 
non-uniform frozen-in field distribution. Also any dynamo mechanism generating 
fields will carry the imprint of inhomogeneous density profile in the star. Thus, 
to maintain the local magneto-static equilibrium, more realistic treatment of the 
equation of state (EoS) of matter for a magnetar requires inclusion of gradually 
increasing magnetic field from surface to center. Massive compact stars are likely 
to develop exotic cores with one possibility of appearance of hyperons with the 
increasing density. In the present work, considering the radial profile of the magnetic 
field and carefully analyzing the different components of the field, we show that the 
pressure of the magnetar matter parallel to the magnetic field exhibits instability. 

2. Model of magnetar matter 

To construct the model of dense matter, we employ non-linear Walecka mean field 
theorjIMI f nU clear matter including the possibility of appearance of hyperons and 
muons at higher density. In the presence of magnetic field, the Lagrangian density 
of the system is C = J2b ^ b + A — jF^F^" , with and Ci are the baryonic 
and leptonic Lagrangian densities respectively in the presence of magnetic fielcP^ 
and the electro-magnetic field tensor. For details, see Ref. [6] 

Total energy density and pressure of the system can be obtained by considering 
the energy-momentum tensor of the system = + T^ u , where and Tf V 
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are the matter and field parts respectively. In the absence of electric field, and in 
the rest frame of fluid 
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when the magnetic field is considered to be along z-direction with B^B^ = —B 2 , B 
is the magnitude of magnetic field, M the magnetization per unit volume, P and e m 
are respectively the pressure and energy density of the matter. This clearly shows 
anisotropic nature of the pressure in the presence of (strong) magnetic field. e m 
is calculated using the charged single particle energy E n — y/p 2 . + m 2 + 2ne\Q\B, 
considering the quantized phase space in the presence of magnetic field, where p z is 
the component of momentum along z-axis, m the mass, e\Q\ the total charge, with 
e being the electron's charge, of the particle, n is the number of occupied Landau 
level. Then matter pressure is P — J2b^b n b + Yli^i n i ~~ £ mi where fi b i and n^; 
are respectively the chemical potentials and number densities for baryons (b) and 
lcptons (I). The density profile of the magnetic field is modeled as^ 
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(2) 



where /3 and 7 are two parameters, rib and no are respectively the number densities 
of matter and nuclear matter, B s and B c are respectively the magnitudes of the 
magnetic field in the surface and center of the underlying magnetar. 



3. Results 

Figure QJi shows EoS for hypernuclear matter in strong and fixed magnetic field 
profiles. For non-zero magnetic field, the pressure splits into the parallel (Py) and 
transverse (Pj_) components, and exhibits anisotropy. It is seen that the low-density 
behavior of EoS with constant magnetic field implies unrealistically large anisotropic 
magnetic field up to the surface of the star, which is inconsistent with the inferred 
surface magnetic field (~ 10 15 G) of magnetars. It is seen from Fig. [TJ) that for a 
given value of /3, the EoS becomes softer with increasing 7. Consequently, beyond a 
certain critical 7 and in a certain density regime, P ceases to increase (and eventually 
decreases) with the further increase in Tib- This implies the onset of instability 
of matter above that value of density for that particular B c and magnetic field 
profile. We also show the results for each j3 with 7 = 1 (minimum value). Note 
that the maximum 7 is taken in such a way that P forms a plateau as a function 
Tib- Furthermore, it is evident that with the decrease of /3, the instability occurs at 
larger values of 7 and n&- 

The instability arises due to the negative contribution from the field energy 
density to the pressure of magnetized baryons and leptons in the direction of the 
magnetic field, which is evident from Eq. ([T]). With the increase of rib, more negative 
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Fig. 1. Left panel (a): Variation of total pressure as a function of normalized baryon number 
density for fixed fields B = (solid line), and B c = 10 18 G and several density profiles: = 10~ 3 , 
7 = 6 (dot-dashed lines), = 10 — 1 , 7 = 4 (dashed lines), and — y co, i.e, B = constant (dotted 
lines). For each pair of curves the upper branch is for P± and the lower branch for Py . Right panel 
(b): Variation of total pressure as a function of normalized baryon number density for different 
magnetic field profiles and B c = 10 18 G. The dots show the reference case B = 0. The solid and 
dashed lines correspond to = 0.1 and 0.001 respectively. For each we choose a pair of 7s; in the 
first case we have 7 = 1 (upper), 4 (lower), whereas in the second case 7=1 (upper), 6 (lower). 

contribution is added to P, and consequently at a certain rib, P ceases to increase 
and then decreases with the increase of rib, rendering instability. 

4. Conclusion 

We have found that for sufficiently large magnetic fields with B c ~ 10 18 G, the mag- 
netar matter becomes unstable. The instability is associated with the anisotropic 
effects arised due to the magnetic field. The onset of instability depends on the 
magnetic field profile and B c , which puts a natural upper bound for the central 
magnetic field of neutron stars, which is 5 x 10 18 G. 

This work was partially supported by the grant ISRO/RES/2/367/f 0-1 f (B.M.) 
and the Alexander von Humboldt Foundation (M.S.). 
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